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VII.  A  statement  of  the  problem  studied  and  a  summary  of  the  most  important  resull 
(The  numbers  refer  to  the  list  of  publications  in  Section  VIII  of  this  report. 
Other  references  are  included  in  the  last  part  of, this  section.) 

The  central  theme  of  this  project  was  the  development  and  analysis  of 
direct  methods  based  on  collocation  for  the  solution  of  singular  integral 
equations  with  a  principal  value  integral.  These  equations  arise  in  such  diverse 
fields  as  linear  elastic  fracture  mechanics,  neutron  transport,  long  water  waves, 
image  reconstruction  and  radiative  transfer.  In  the  classical  approach,  the 
singular  integral  equation  is  regularized  to  yield  a  Fredholm  integral  equation 
of  the  second  kind.  The  numerical  implementation  of  the  regularization  is  usually 
quite  cumbersome.  While  direct  methods  proposed  by  Erdogan(1969) ,  Erdogan 
and  Gupta(1972)  ,  and  Theocaris  and  Ioakimidi s ( 1 977 )  have  been  used  in  Fracture 
Mechanics,  their  convergence  and  stability  has  received  scant  attention.  The 
work  at  Stony  Brook  has  made  important  strides  in  both  directions. 
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We  have  analysed  the  behavior  of  approximate  solutions  enhancing  our  under¬ 
standing  of  the  convergence  and  stability  of  the  methods  based  on  Gaussian 
quadrature.  We  have  also  proposed  methods  which  can  be  employed  in  the  situations 
where  the  Gaussian  quadrature-collocation  schemes  fail. 


The  key  to  successful  analysis  of  the  aforesaid  direct  methods  was  the 
discovery  that  the  coefficient  matrix  of  both  Gauss-Chebyshev  and  Lobatto- 
Chebyshev  techniques  have  a  closed  form  inverse,  which  is  the  transpose  of  the 
coefficient  matrix  multiplied  by  a  diagonal  matrix  [1],  [2].  Using  Jackson's 
theorem  on  "best  approximation",  the  convergence  of  interpolator  polynomials 
can  be  easily  established.  By  refining  the  techniques  we  are  able  to  prove  [3] 
the  convergence  uuder  weaker  condition.  [4]  contains  an  extension  of  the  analysis 
of  [l]  to  the  complete  singular  integral  equation.  The  nodes  of  quadrature  and 
collocation  now  are  the  zeros  of  Jacobi  polynomials.  To  cover  the  general  case 
of  a  Fredholm  kernel  added  to  a  Cauchy  kernel,  the  theory  of  collectively  compact 
operators  is  used.  [5]  is  addressed  to  the  problem  of  computing  bounded  solutions. 

In  this  case  a  compatibility  condition  has  to  be  satisfied.  Often  the  discrete 
system  representing  the  continuous  is  overdetermined  and  has  one  equation  more 
than  the  number  of  unknowns.  We  find  that  the  solutions  vary  with  the  equation 
ignored  (the  usual  engineering  practice)  and  the  error  peaks  in  the  vicinity  of 
the  collocation  point  which  is  discarded  to  obtain  a  compatible  system.  Use  of 
general  collocation  points  in  conjunction  with  Gaussian  quadrature  is  examined 
in  [6].  The  question  of  an  optimal  selection  of  nodes  is  still  open. 

While  both  Gauss-Chebyshev  and  Lobatto-Chebyshev  methods  are  easy  to  code, 
they  are  not  suitable  for  functions  which  are  highly  oscillatory  in  some  small 
neighborhood.  Since  the  collocation  points  are  the  zeros  of  Chebyshev  polynomials, 
the  only  way  their  concentration  can  be  increased  is  by  increasing  the  degree  of 
the  polynomial.  Since  the  largest  distance  between  the  zeros  of  Chebyshev  polynomials 
is  approximately  ir/n,  for  a  resolution  of  the  order  .001  a  polynomial  of  order 
3000  has  to  be  considered. 

In  [7]  we  propose  a  method  based  on  the  use  of  piecewise  linear  functions. 

[8]  develops  the  procedure  further  employing  cubic  splines.  The  order  of  singularity 
is  explicitly  built  into  the  solution.  Na.iet-Nasser  and  Horii  (1982)  found  that 
the  results  obtained  through  piecewise  polynomial  approximation  indistinguishable 
from  the  corresponding  results  obtained  by  the  Erdogan-Gupta  method  for  the  problem 
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of  kinked  crack  extension  under  tension.  [9]  contains  a  review  of  the  Gaussian 
integration  methods  for  the  numerical  solution  of  singular  integral  equations. 

m 

Some  of  the  results  of  our  investigation  are  yet  to  be  published.  To  reduce 
an  integral  equation  to  an  algebraic  system  whose  solution  leads  to  a  polynomial 
approximation  to  the  solution  usually  one  of  the  two  approaches  is  adopted:!.  The 
weighted  norm  of  the  difference  between  the  left  hand  side  and  the  right  hand  side 
is  minimized  by  choosing  the  coefficients  of  the  polynomial  (The  Galerkin  Method) 
2.  The  equation  is  satisfied  only  over  a  finite  set  of  points  and  the  integral 
is  replaced  by  a  Quadrature  formula  (The  collocation  method).  We  have  been  able 
to  show  that  the  two  solutions  are  related  by  means  of  a  Gaussian  quadrature 
formula  [10]  .  [11]  contains  a  comparison  of  the  accuracy  of  stress  intensity 
factors  for  a  Griffith  crack  problem  using  Gauss-Chebyshev  and  Lobatto-Chebyshev 
quadrature  formulas.  [12]  contains  a  proof  that  the  coefficient  matrix  obtained 
through  the  "tanh  rule"  and  collocation  is  nonsingular.  Finally,  [13]  treats 
the  problem  of  a  cruciform  crack  opened  by  internal  pressure  varying  with  time. 

For  time  dependent  problems  it  brings  into  focus  the  need  for  methods  which  are 
both  fast  and  accurate. 

We  are  in  the  process  of  writing  the  final  versions  of  the  papers  [10,11,12] 
and  U.S.  Army  Research  Office  will  be  furnished  copies  when  they  are  ready. 
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